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‘Preferential concentration’ {Phys. Fluids A3:1169-78, 1991) refers to the clustering of 
inertial particles in the high-strain, low-rotation regions of turbulence. The ‘centrifuge 
mechanism’ of Maxey (J. Fluid Mech. 174:441-65, 1987) appears to explain this phe¬ 
nomenon. In a recent paper, Bragg & Collins {New J. Phys. 16:055013, 2014) showed 
that the centrifuge mechanism is dominant only in the regime St 1., where St is the 
Stokes number based on the Kolmogorov time scale. Outside this regime, the centrifuge 
mechanism gives way to a non-local, path-history symmetry breaking mechanism. How¬ 
ever, despite the change in the clustering mechanism, the instantaneous particle positions 
continue to correlate with high-strain, low-rotation regions of the turbulence. In this pa¬ 
per, we analyze the exact equation governing the radial distribution function and show 
how the non-local clustering mechanism is influenced by, but not dependent upon, the 
preferential sampling of the fluid velocity gradient tensor along the particle path-histories 
in such a way as to generate a bias for clustering in high-strain regions of the turbulence. 
We also show how the non-local mechanism still generates clustering, but without pref¬ 
erential concentration, in the limit where the timescales of the fluid velocity gradient 
tensor measured along the inertial particle trajectories approaches zero (such as white- 
noise flows or for particles in turbulence settling under strong gravity). Finally, we use 
data from a direct numerical simulation of inertial particles suspended in Navier-Stokes 
turbulence to validate the arguments presented in this study. 

Key words: 


1. Introduction 


Inertial particles that are initially uniformly distributed throughout an incompressible 
turbulent flow will develop dynamically evolving spatial clusters. Early studies observed 
that inertial particles not only cluster but preferentially concentrate in strain-dominated 


regions of the turbulent fluid velocity field, away from regions of strong rotation (Maxey & 

Corrsin|1986| Maxey|19871 

Squires & Eaton|1991[|Eaton & Fessler|1994||Wang & Maxey 

1993 

Sundaram & Collins 

1997 

). Here and throughout we make a distinction between 


the terms ‘clustering’ and ‘preferential concentration,’ where the former refers to the non- 
uniform distribution of the particles in space, irrespective of where they cluster, and the 
latter refers to the clustering of particles in high-strain, low-rotation regions of the flow. 
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The reason for making this distinction is that clustering can occur without preferential 
concentration, a scenario that we discuss later in the paper. 


The traditional explanation for the mechanism of clustering is due to Maxey (1987), 


who argued that particles are centrifuged out of regions where the fluid streamlines 
exhibit strong curvature into regions of high-strain. This physical argument, aside from 
being intuitive, explains the preferential concentration effect. However, over time the 
completeness of this explanation for clustering has been questioned. For example, |Bec| 
(2003) described simulations of particles in white-in-time random flows that exhibited 


strong clustering, yet the centrifuge mechanism does not operate in such flows. 
In a recent paper. 


Bragg & Collins (2014) considered in detail the physical mechanism 


responsible for the clustering of inertial particles in the dissipation range of isotropic 
turbulence. We showed that in the regime St 1 (where St = Tplr^ is the Stokes number, 
Tp is the particle response time and Tp is the Kolmogorov timescale), the mechanism for 


clustering predicted by the Zaichik & Alipchenkov theory (Zaichik & Alipchenkov 


2007 |2009 hereafter ZT) is the same as that predicted by the Chun et al. theory (Chun 


2003 


et aL||2005[ hereafter CT), which is essentially an extension of the classical centrifuge 
mechanism of Maxey ( |1987 ). As St is increased beyond this regime (i.e., St ^ 0(1)), 
the system undergoes a bifurcation and the particle velocity dynamics become non-local 
(by ‘non-local’ we mean that the particle velocity at a given time depends upon its 
interactions with the turbulent fluid velocity field at earlier times along its path-history), 
and we showed that this non-locality contributes to the inward drift mechanism described 
by the ZT through a path-history symmetry breaking effect. Figure [T] illustrates two 
particle pairs passing through separation r at time t. For St ^ 0(1), the particles retain 
a finite memory of their interaction with the turbulence along their path histories, over 
a timeframe 0{t — Tp) ^ s Si t. The particle pair separation is given by rP{s\r,t), where 
rP is the time dependent separation, and \r,t denotes that the pair has separation r 
at time t. For the solid-line trajectory in figure rP{s\r,t) ^ r and for the dashed-line 
trajectory, rP{s\r, t) ^ r. Since the fluid velocity difference. Ait, on average increases with 
separation then particle pairs with rP(s|r, t) ^ r will have experienced, on average, larger 
values of Am in their path-history as compared with particle pairs with rP{s\r,t) ^ r. 
This asymmetry in the particle pairs’ path histories breaks the symmetry of their inward 
and outward motions and gives rise to a net inward drift, leading to clustering. 


We also argued that the theory for clustering in Gustavsson & Mehlig (20111 is qual¬ 


itatively similar to that in the ZT in the regimes St 1 and St ^ 1. The non-local, 
path-history symmetry breaking mechanism in the ZT and the non-local ‘ergodic mech- 


anism’ in 

Gustavsson & Mehlig 

(2011 

) both explain the clustering in white-in-time flows 

observed 

in 

Bee ( 

2003 

). We will return to this point later in the paper. Note that the ZT 


is more general in its description, as it applies to all St and correlation timescales for the 


fluid, whereas Gustavsson & Mehlig (2011) treat these timescales as a small parameter. 


There is however an outstanding issue that remains to be explained. |Bragg fc Collins| 
(2014) showed that the non-local clustering mechanism in the ZT makes the dominant 
contribution when St = 0(1) (see figure 3 in Bragg fc Collins|2014 and the surrounding 


discussion). However, we also know from numerical simulations (e.g.. Squires & Eaton 
[1991 ) that at 5't = 0(1) there exist strong correlations between the location of the clusters 


and the local topology of the fluid velocity gradient field, i.e., preferential concentration. 
In the regime S't <C 1, the phenomenon of preferential concentration is easily explained 
since the mechanism generating the clustering depends upon the preferential sampling of 
the local fluid velocity gradient ( Ghun et a7|2005 Balachandar fc Eaton|2010 ). However, 
for St = 0(1) the clustering mechanism is strongly non-local, and this symmetry-breaking 
clustering mechanism does not explicitly depend upon the local flow topology. How then 
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Figure 1: Illustration of two particle paths arriving at separation r at time t. The solid 
lines represent a pair of particle trajectories with separation r^’(s|r, t) ^ r and the dashed 
lines a pair with rP{s\r,t) ^ r. 

does the non-local clustering mechanism continue to generate preferential concentration? 
In this paper, we address this question through an analysis of the exact equation for the 
radial distribution function (RDF), an important measure of clustering. In particular, 
we resolve the issue by showing that the non-local clustering mechanism is influenced 
by the preferential sampling of F, the fluid velocity gradient tensor, along the path 
histories of the particles in such a way as to generate stronger inward drift velocities into 
regions where F at the current particle position exhibits high-strain and low-rotation. Our 
analysis also shows how in flows where the centrifuge mechanism does not operate and 
the preferential sampling effect is absent, the non-local mechanism generates clustering 
without preferential concentration. 


2. Theoretical Analysis 

2.1. Derivation of the drift velocity conditioned on the local flow structure 

We consider the relative motion between two identical point particles, a ‘primary’ particle 
and a ‘satellite’ particle, the primary particle being the reference particle. The particles 
experience Stokes drag forces, but do not interact with each other through physical col¬ 
lisions or hydrodynamic interactions and are assumed to be at low enough concentration 
to not affect the turbulence (i.e., ‘one-way coupling’). The turbulence is assumed to 
be statistically stationary, homogeneous and isotropic, and gravitational settling is ne¬ 
glected. The particle-to-fluid density ratio is large, corresponding to a gas-solid system. 
Under these assumptions, the equation governing the relative motion of the two particles 
reduces to ( Maxey fc Riley |[l983| ) 



( 2 . 1 ) 
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where (t), {t), (t) are the particle pair relative separation, velocity and accelera¬ 

tion vectors, respectively, xP{t) denotes the primary particle position and Au{xP{t), rP{t), t) 
is the difference in the fluid velocity evaluated at the positions of the two particles. Since 
we are interested in the case where u{x, t) is statistically homogeneous, we shall drop 
the x^it) argument when Au appears in statistical expressions. 

For the system governed by (13 the exact equation governing the probability density 
function (PDF), 

p(r, w, t) = — r)S{wP{t) — w)'^, 

that describes the distribution of rP{t),wP{t) in the phase-space r,w is 


dtP = -Vr-pw + {StT,^) ^ pw - (StT^,) ^ p(^Au{rP{t),t)'j 


( 2 . 2 ) 


where {■)r.w denotes an ensemble average conditioned on T'^’(t) = r and wP{t) = w. A 
commonly used statistical measure of particle clustering is the radial distribution func¬ 
tion (RDF, McQuarrie|1976 1, which is defined as the ratio of the number of particle pairs 
at separation r to the number that would be expected if the particles were uniformly 
distributed. An exact equation for the statistically stationary RDF, g{r), can be con¬ 
structed by multiplying the stationary form of (2.2) by w and then integrating over all 
w yielding 


0 = g(^Au{rP{t),t)'j^ - Sir, 


qp 

r ]^2 


Vrg- StTriP Vr- S^, 


where 


gir) = 


N{N-1) 


p{r, w) dw, 


(2.3) 


(2.4) 


N is the total number of particles lying within the control volume V, n = N/V is 
the number density of particles, and ^'^(t’) = {w^{t)wP{t))r is the second-order particle 
velocity structure function. 

Equation ( |2.3[ ) is exact, but unclosed, due to the term {Au{rP{t),t))r which describes 
the average of Am experienced by particle pairs at the separation T'^’(t) = r. Even though 
{Au{r,t)) = 0 for isotropic flows, {Au{rP{t),t))r ^ 0 for finite St because the inertial 
particles preferentially sample the underlying field Au{x,r,t). In the ZT, this term is 
closed by approximating Au{x,r,t) as a spatio-temporally correlated Gaussian field and 
using the Furutsu-Novikov closure method ( Zaichik &: Alipchenkov|2007 [Bragg fc Collins 
2014). However, for the purposes of this analysis, it is not necessary to introduce any such 


closure approximations. Recall that the purpose of this analysis is to demonstrate that 
the non-local clustering mechanism, which dominates at St = 0(1), is consistent with 
the phenomenon of preferential concentration. It is unnecessary to close {Au{rP(t), t))^ 
for two reasons. First, Bragg & Collins ( 2014| ) showed that at St = 0(1) the drift velocity 
contributions coming from {Au{rP{t),t))r are negligible compared to the drift mecha¬ 


nism associated with StTrj Vt 


S '2 in (2.3). The reason for this is that at St = 0(1) in 


the dissipation range of the turbulence, the particle relative velocities exhibit ‘caustics’ 
(Wilkinson & Mehlig 2005 Bee et al. 2010 Salazar & Collins 20126), which implies 


they are much larger than the fluid relative velocities. A consequence of this is that 
contributions to the particle relative motion coming from their local interaction with 
the turbulence (e.g. Au{rP{t),t)) become negligible. Second, any contributions to the 
clustering mechanism associated with {Au{rP{t),t))r. are manifestly consistent with the 
phenomenon of preferential concentration since this term arises because of preferential 
sampling effects. 



























Non-local clustering mechanism and preferential concentration 


The challenge therefore is to show that the drift velocity Vt 


5^2) which generates 


the clustering and is strongly non-local at St = 0(1), is consistent with the phenomenon 
of preferential concentration. The physical interpretation of StTn Vr* has already 
been discussed in the introduction. We refer the reader to Bragg & Collins (2014) for a 
more detailed explanation. 

Since we are interested in clustering in the dissipation range we take 


Au{xP{t),rP{t),t) = T{xP{t),t) • rP{t), 


where r(a;, t) = 'Vxu{x, t) is the fluid velocity gradient tensor. If the non-local clustering 
mechanism for particles with St = 0(1) is consistent with the phenomenon of preferential 
concentration, then it must be the case that contributions to Str^^ V^.- are larger for 
particles that are (at present time, t) sitting in high-strain, low-rotation regions than 
they are for particles sitting in high-rotation, low-strain regions of the flow field. 

In order to formally consider this possibility, we expand the phase-space description 
of the drift velocity as shown below 


- StTr^ Vr-SP = 



Drift velocity conditioned on t) — Z 


dZ, 


(2.5) 


where Z is an independent phase-space variable, q{Z) = {S{ZP{t,t) — Z)), 

ZP{t',t'') = S{xP{t'),t') : S{xP{t''),t'') —'J^{xP{t'),t') : 7?.(a;*’(t"), t"), S and Tt are the 
fluid strain-rate and rotation-rate tensors, respectively, defined as 5 = (l/2)(r -f- F^), 
TZ = (l/2)(r — r^), and {■)r,z denotes an ensemble average conditioned on rP{t) = r and 
ZP{t,t) = Z. We have defined ZP in terms of two time arguments t' and t" since quan¬ 
tities appearing later in the analysis will involve the more general expression ZP{t',t''). 


Through (2.51 we are able to consider how the non-local contribution to the inward drift 


velocity is related to the local properties of F at the particles current location xP{t), 
where the relevant properties of T{xP{t),t) are characterized by its invariant ZP{t,t). 

We shall now define what we mean by high-strain, low-rotation regions and high- 
rotation, low-strain regions. A rigorous definition of these regions could be given in terms 


of the invariants and eigenvalues of T{x,t) (e.g. Chong et al. 1990 Rouson & Eaton 
|2001[ Salazar & Collins 2012a), however for our purposes this is not necessary. In the 
regime St <C 1, the single particle velocity vP{t) can be approximated as being a field, i.e. 
vP{t) « v{xP{t),t) with a divergence ^x‘v{xP{t),t) « —StTr^ZP[t,t) implying particles 
will cluster in regions where ZP{t,t) > 0 ( Maxey|[l987 ), which are commonly referred to 
as high-strain, low-rotation regions. Numerical simulations have not only confirmed this 
for St 1, but have shown that inertial particles continue to preferentially concentrate 
in ZP{t,t) > 0 regions even when St — 0(1). Since it is the clustering behavior in these 
regions that is of interest, what we want to demonstrate is that for St = 0(1) 


— StTrj Vr 


wP(t)wP{t)\ < —StTrj Vr- (wP {t)wP {t)\ 

/ \ / r 


.zi-i 


( 2 . 6 ) 


StTjj Vr 


+ 00 

wP(t)wP(t))^^^^^= f q{Z)\ - StTr^V, 


wP (t)wP (t) 


dZ, 


where 
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and 


—Str^ Vj 


wP{t)wP{t) 




^ / 0(2)1 

-5tT,V,.(M;P(t)M;P(t)) 

J 1 

— OO 

\ / r,Z 1 


dZ. 


Note that by definition 

— StTn 'Vr‘ So = —StT„ 'Vr‘ 


wP{t)wP(t) ) - StTrf Vr- 


wP{t)wP{t) 


r,Zp 


The formal solution for wP{t) for particles in the dissipation range governed by (2.1) 
is (ignoring initial conditions since we are interested in the statistically stationary state) 


wP{t) = {StT^)-^ j • rP{t')dt', 

0 

and using this we obtain 


(2.7) 


— StTri Vr 


wP {t)wP (t) 


= -{StTrjY 


V^- Qdt'dt", (2.8) 


0 0 


where the components of Q are 


Q,, (r, Z, t, t', t") = {TU=^P(t'), t')rP^{t')V,r.{xP{t"),t")rP^{t") 


r.Z 


(2.9) 


We now seek a simplification of (2.9) that is valid in the regime of interest, namely 
St = 0(1) and r t], where rj is the Kolmogorov length scale. Since we are interested 
in St = 0(1), neither St nor 1/St can be treated as a small parameter. However in this 
regime, the particle velocity dynamics are strongly non-local, implying wP 3> Am, which 


is associated with the formation of caustics in the particle velocities (Wilkinson & Mehlig 
2005 Bee et al.|2010[ Salazar &: Collinsj2012& ). This suggests that we can define the small 


parameter ^ = \{Au{r,t)Au{r, t))|/|S' 2 (T’)|, which in the caustic regions satisfies fi 1. 
We therefore introduce the perturbation expansion 

Q = Q[0l+^g[ll+^2g[2]^,,.^ 


where the components of are 


( 2 . 10 ) 

( 2 . 11 ) 

/ r \ ! ^ 

This expression for follows from the fact that in the limit /r —)■ 0 the particle motion 
approaches ballistic motion, implying that the correlation between rP and F vanishes, and 
the particles uniformly sample F. In the regime /j, ^ 1, to leading order Q « + 

which can be used in (2.8) to construct an expression for —StTjj {wP{t)wP{t))r,z- 
The term {rP{t')rP{t"))r in (2.11) is related to the average growth of the particle pair 
separation rP backward in time (noting that t' ^ t and t” ^ t), evaluated along pair 
trajectories that satisfy rP{t) — r. For t' = t" the quantity is the backward-in-time mean 
square separation of the particle pair. The quantity may be expressed as 

(rP{t')'i'P{t")^ = rr + &, (2-12) 

where @{r,t,t',t”) is the backward-in-time separation tensor, satisfying @{r,t,t,t) = 0 . 
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In the regime /i <?; 1, using the formal solution for for particles governed by (13 , we 
obtain the approximation (see Bragg et aL||2014 1 


0 = {StT^fSl 


I _ g(*-*')/'fp _ g(*-*")/'rp _j_ g(2i-t'-t")/Tp 


+ 0{p). 


(2.13) 


The result in (2.13) is only valid for max[t — f,t — f] ^ 0{Tp) since for larger time 
separations the process transitions into a different regime ( [Bragg et a7||2014 ). However, 
particles only retain a memory of their interaction with the turbulence up to times 


0{t — Tp) in their path-history, making the approximation in (2.13) sufficient to capture 
the dominant contribution to the drift velocity. 


Substituting (2.12) into (2.11), and taking the divergence yields 


= Vr^\rmrn]\^i7n{x^{t'):t')Tjn{x^{t"),t") 

+ Vp^[0„„](r,^(a;P(t'),t')r,„(a:^(O,O)^+O(M). 


(2.14) 


Using (2.13) in this expression, and substituting this into ( |2.8[ ) and simplifying yields 

t t 

StTpVr- (wP[t)wP{t)^ ^=-{StTp)-^^ f 


0 0 




0 0 


dt'dt" 


(2.15) 

where yP{f,t") = S{xP{t'),t') : S{xP{t''),t") -^n{xP{f),t') ■.n{xP{t"),t"), 

0|| (r, t, t', t") = r~^rr : 0(r, —t, —t', —t”), and we have introduced the approximation in 
0 that Sf, « <S' 2 ||I, which is valid for St > 0(1) (Wang et aL|2000). The notation Z[+’“l 


in (2.15) denotes that the conditionality is either to be taken as Zl+l or ZH. For an 


isotropic system, the drift velocity between the particles is given by the projection of 
(2.15) along r~^r such that the drift velocity only depends on r = |r| and not the vector 


r. Consequently, the sign of r in (2.15) is irrelevant in the following discussion. 

The physical interpretation of terms ® and @ may be explained as follows. Suppose 
that during the particle timescale, Tp, the pair separation remains almost constant. 
In this case term @ would vanish (because 0|| and Vp0|| are zero if rP is constant, see 


(2.12)) but term 0 would be finite and describes a drift velocity arising because the 
primary particle preferentially samples F along its trajectory. This scenario occurs in the 
regime St 1, as Tp is so small that rP remains nearly constant over the timescale of the 
particles. Consequently, term 0 ^ term (0 and the inward drift arises from the biased 
sampling of F along the path-history of the particle for times up to 0(t — r^). In this 


regime, the sum of term 0 for Zl+l and 1 yields the CT result (Chun et al. 
Bragg fc ColIh^lMIt 


2005 


-StTp Vr- SP = -StTj-{^ZP{t,t)). 


For St = 0(1), in contrast, the pair separation rP changes significantly over the timescale 
of the particles. This change in their relative separation affects the fluid velocity differ- 
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ences the particles experience along their path-history since Au{xP (t), (t), t) oc rP{t) in 

the dissipation range, and this gives rise to the path-history symmetry breaking effect de¬ 
scribed earlier. The time depen dent r elative separation between the primary and satellite 
particles is described by 0 (see (2.12)), which appears in term (2) through ©y. The growth 
of 0 backward in time ensures that particle pairs arriving at r, t from larger separations 
dominate the behavior of {w'P[t)w'P{t))^ relative to particle pairs arriving at r,t 

from smaller separations. An important point is that although St = 0(1) and Am oc r 
gives rise to the operation of the non-local, path-history symmetry breaking mechanism, 
the strength of the resulting inward drift in the dissipation range is influenced by the 
way the particles have interacted with S and 'R. along their path-history. In term 
this dependence is described by {y^{t',t")) 2 .[+,-\- We will show that it is the fact that 
inertial particles interact with the fields S and TZ differently that allows the non-local 
clustering mechanism to generate stronger inward drift velocities into high-strain regions 
of the turbulence. 


2.2. Analysis of the drift velocity for finite flow time scales 
Let us now consider how preferential sampling of F along the path-history of the particles 
influences terms ( 1 ) and and in particular, how this affects the strength of the drift into 
ZW and regions. In what follows, we assume that the autocorrelations associated 
with S{xP{t),t) and 'TZ{xP{t),t) are positive (at least for max.[t — t',t — t"] < 0(rp)). 


This is known to be true for St = 0 particles (e.g., Girimaji & Pope 

1990 

) and it seems 

reasonable to assume that this is also true for St > 0 particles (in |3 

we will show using 


DNS data that this is indeed the case). 


Term 0 in (2.15) contains 


zp{t',n)^^^_^ = (^s{xP{fl),t'):s{xP{nx)-n{xPifl),t'):n{xP{n,n) 


Zl+'-l 


For ZA] the contribution involving S dominates, and for Z^ 1 the contribution involving 
IZ dominates. Consequently 




(2.16o) 

(2.166) 


which leads to the expectation that the mechanism associated with term 0 causes the 
particles to drift out of ZA] regions and into Zl+1 regions. Therefore term 0 is consistent 
with preferential concentration. 

An implicit assumption made in ( 2.16[ ) is that the expression {Z'P{t',t'')) 2 : 1 +,-] remains 
of the same sign for max[t — t',t — t"] ^ Oirp). That is, if a particle is in a region 
at time t, then it has, on average, been in a Z^^ region for as long as it can remem¬ 
ber its interaction with the turbulence along its path-history, and similarly for ZA], 
As an estimate, this is satisfied provided that the Lagrangian timescales of S{xP{t),t) 
and TZ{xP{t),t), ts and r-ji, respectively, satisfy ts 0{Tp) and t-ji 0{Tp). DNS data 
( Ireland et al.|[2M4 ) shows that this condition is satisfied for St < 0(10). Another, asso¬ 
ciated reason why the conditions ts ^ 0(Tp) and t-ji 0(Tp) must be satisfied is that if 
Ts <C Tp and T-ji <C Tp then the dominant contribution to the drift velocity would be from 
the particles’ interaction with F in regions that are uncorrelated with F at its current 
position. In this situation the bias in the drift velocity for given values of Z would nec¬ 
essarily vanish, as would the preferential concentration, a scenario that we shall return 
to later. 














Non-local clustering mechanism and preferential concentration 
Term (2) contains 

y^t',t "))= (s{x^{t'),t') : S{x^{t"),t") + n{x^{t'),t') : n{x^{t”)X)) 


Similar to {Z^{t',t")) 2 .[+-], the contribution involving S dominates i")) 2 [+i, and 

the contribution involving 7L dominates {y^{t',t"))z[-]- Unlike {ZP{t',t"))z[+-] how¬ 
ever, {ypft',t")) 2 ,[+,-] is of the same sign for and and therefore the non-local 
clustering mechanism associated with term (2) causes the particles to drift into both ZUI 
and regions. Nevertheless, since inertial particles with St = 0(1) preferentially sam¬ 
ple strain-dominated regions of the flow and under sample rotation-dominated regions 
(e.g., Salazar fc Collins||2012a) then we expect that 




(2.17) 


Consequently, based on the predicted inequalities in (2.161 and (2.17), and since VrOii > 0, 


(2.15) leads to the expectation that even for St = 0(l) when the clustering mechanism is 
strongly non-local, the particles continue to drift more strongly into regions, consis¬ 
tent with the phenomenon of preferential concentration. The physical explanation is that 
because the inertial particles preferentially sample F along their trajectories, the fluid 
velocity differences driving the particles together are experienced by the particles to be 
greater in regions than in Z^~^ regions. This preferential sampling is a consequence 
of the particles being centrifuged away from regions of strong rotation along their tra¬ 
jectories. Nevertheless for St = 0(1), the centrifuge mechanism is not the primary cause 
of clustering, but rather it simply influences the clustering behavior through the way it 
causes particles to preferentially sample F along their path-histories. A consequence of 
this is that St = 0(1) particles can cluster even in situations where the centrifuge mech¬ 
anism does not operate, for which biased sampling and hence preferential concentration 
disappear. We now consider this situation. 


2.3. Analysis of the drift velocity for white-noise flows 

In the limit where —>■ 0 and —>■ 0, such as in white-in-time random flows or for parti¬ 
cles falling through turbulence in the limit of strong gravity ( Bee et qZ.| 2014 Gustavsson 
fc Mehli^|2014 |, the centrifuge mechanism does not operate and there is no preferential 
sampling of F. In this limit 


ZP(t',t") 


ZP{t,t) 


Zi+'-'i 


/zl+.-l’ 

t" = t' = t, 

-n, 

t" <t,t' < 


where 


zifl)= (zp it', n) 

J —oo ' 


dt" 


(2.18) 


(2.19) 


The t' < t, t" < t result follows from the fact that in the limit 0 and t-ji —>■ 0, 

(ZP(t',t"))zi+-] = 0 if ^ t", and since ZP(t',t') and ZP(t,t) would be uncorrelated in 
this limit then {ZP(t' = (ZP{t',t')) for t' ^ t. Furthermore, in the statistically 

stationary state, {ZP(t',t')) = (ZP(t,t)) = 0 Z{t') = 0 and 


ZP(t,t) 


zw 
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such that the contribution from term (1) to the total drift velocity is zero. However, in 
the same limit 


ynt\n 


ynt,t) 




) zi+’-S 

t” = t' = t, 

-n, 

t” <t,t'< t. 


where 


J —oo ' 


dC 


( 2 . 20 ) 


( 2 . 21 ) 


Unlike {ZP{t', t')), {yP{t',f)) ^ 0 ^(^0 7 ^ 0 in the statistically stationary state. The 

result in (2.20) for t" = t' = t does not render term @ dependent upon Z in the limit 
—)• 0 and —)■ 0 because in the integrand of term (y^ (t', t")) ^ [+, _] is m ultiplied 

by Vr0||, and for t” = t' = t, VrOy = 0. Consequently, since the result in (2.20) is inde¬ 
pendent of Z for t" <t, t' < t, then in the limit —)■ 0 and t-ji —>■ 0, term (3) generates 
a finite contribution to —StTri S 2 but without a bias towards any particular region 
in Z space. That is, term (2) becomes independent of Z and hence the drift velocity is 
equally strong for all Z. 

2.4. Comparison with Custavsson & Mehlig (2011) 

The result presented in §2.3| is consistent with the results derived by Custavsson & 
Mehlig (2011), who showed that in the limit —>■ 0 and Tn —7 0, the particles clus¬ 


ter, but without the preferential concentration effect. As discussed in [Bragg fc Collins 
(|2014|), the ‘non-ergodic’ and ‘ergodic’ clustering mechanisms described in Custavsson 
fc Mehligj ( [Ml] ) are qualitatively similar to the clustering mechanisms described by 
—StTrj Vr- S 2 in the regime Tp min[T 5 ,T 7 ^] and Tp 3> max[T 5 ,T 7 j], respectively. The 
terms non-ergodic/ergodic in this context refer to the presence/absence of preferential 
sampling effects. Because of its perturbative construction, the clustering theory in |Cus-| 
tavsson & Mehlig (2011) was not able to describe the interplay between the non-ergodic 
and ergodic clustering mechanisms for St = 0(1) and Ku ^ 1 (where Ku is the non- 
dimensional Kubo number, which is ~ 1 for real turbulence, see Duncan et al. 2005), 


although their simulations indicate that both make a substantial contribution to the 
clustering in this regime. The theoretical analysis in this paper demonstrates how these 
two effects contribute in this regime: The non-local, path-history symmetry breaking 
mechanism, which operates even in the purely ergodic limit, is nevertheless influenced 
by non-ergodic effects when St = 0(1) (and Ku ~ 1), and these non-ergodic effects are 
what allow the non-local clustering mechanism to generate preferential concentration in 
this regime. 

2.5. Validity of dissipation range scaling 

In closing this section we briefly comment on the validity of assuming dissipation range 
scaling for Au{x,r,t). Since depends non-locally in time on Am, then by introduc¬ 
ing Au{xP(t), rP{t), t) = T{x^{t),t) • r^ft) into ( | 2 . 1 | ) we are implicitly assuming in our 
analysis that this assumption is valid up to times Oft — Tp) along the path-history of 
the particle pair. The dissipation range scaling Au{x,r,t) oc r is supposed to be only 
valid for r <C ?? ( Pope|20M ). This could be problematic for our analysis of particles with 
St = 0(1) that may retain a memory of their interaction with the turbulence at \rP\ > 77 . 
However, the problem is mitigated by the fact that numerical studies of Navier-Stokes 
turbulence have shown that Au{x, r,t) oc r up to r = O{10ri) (Ishihara et al. 2009). Thus, 
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the effective dissipation length scale is an order of magnitude larger than the Kolmogorov 
length scale, likely due to nonlinear depletion effects in the turbulence ( Frisch|1995 1 that 
allow the viscosity to act on scales an order of magnitude larger than would be expected 
based on a scaling argument. St = 0(1) particles at rj have a temporal memory 

that is sufficiently small that interactions with the turbulence at separations greater than 
0 ( 1077 ) in their path-history only weakly affect the drift, an argument supported by the 


results in Ray & Collins (2013). 


3. Tests using DNS 

In this section, we use DNS data to test the arguments presented in the previous 
section. In particular, we will calculate from DNS data 


ZP{t',t") 


zi+’-l 


ynt',n 


Zl+'-l 


along with S’^n (required in 0||), evaluate the integrals in (2.15), and then compare the 


strength of the drift velocity in the and 1 regions. The DNS simulations were of 
statistically stationary, homogeneous and isotropic turbulence at Re\ = 224, where Re\ 


is the Taylor microscale Reynolds number. Details on the DNS can be found in Ireland 


et al. (2013). 


In computing the results, we replace the lower integral limit of 0 in (2.15) with t — Tp. 


The reason for doing this is that the closure approximation used for 0 is only valid 
for max[< — t', t — t"] < 0(rp), and gives unphysical results at larger time separations. 
However, the cut-off is sufficient to capture the dominant contribution to the drift velocity 
since, as discussed in the particles only retain a memory of their interaction with the 
turbulence for times up to 0{t — Tp) along their path-history. 

To make the notation in what follows more succinct we define 

= —StTp Vr’ (w^ {t)wP (t)'^ 

and since we are considering isotropic turbulence it is only the parallel component of 
the drift velocity that is of interest, = r~^ r- We also use the additional 

subscript of either 1 or 2 (e.g. o[|~^J) to denote that we are considering the contribution 




to olr*"’ ^ coming from either term (I) or term ( 2 ) in (2.15). 


In order to examine whether the drift velocity contains a bias for or zi-1 we 
consider 


r _ ^[-] /^[+] 

T _ .[-] /^[+] 

J 2 = 0j|,2/®||,2- 


(3.1) 

(3.2) 


Since the r dependence of the drift velocity contributions to terms ® and (2) are the same 
for and Z^~\ R and J 2 are functions of St, but not of r. As discussed earlier, the 
sum of the 2^1+1 and Z^~^ contributions from term (I) is zero in the absence of preferential 
concentration, implying that in the absence of preferential concentration R = —1. For 
term @ however, the absence of preferential concentration is denoted by R = 1, since, 
as demonstrated in (2.2 in the absence of preferential concentration, term @ becomes 
independent of Z. The results for Ji and R are shown in hgure[^ For St > 0.7, R and 
R are evaluated using (12.151); for St ^ 0.2, we also evaluate R and R using the St 1 
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Figure 2: Plots of (a) Ji and (b) J 2 as a function of St. 


approximations 


J +,-1 




()'-l - c)t+l 


0 


(3.3) 

(3.4) 


(such that J 2 is undefined and is therefore not plotted for St < 0.2). The results in figurej^ 
show that for St > 0, Ji > —1 and J 2 < 1, consistent with preferential concentration. As 
St is increased from 0, Ji initially increases, reaches a maximum at St « 0.7 and then 
begins to decrease again, consistent with known DNS data for the preferential concentra¬ 
tion effect (e.g., see hgure 1 in Bee et fl/)]|2007 |. That — 1 ^ Ji <0 reflects the fact that 
the clustering mechanism associated with term (1) causes particles to drift away from 
regions and drift into and is directly associated with the centrifuge mechanism act¬ 
ing along the particle trajectories. In contrast, 0 ^ J 2 < 1 signifying that the mechanism 
associated with term the non-local clustering mechanism, causes particles to drift into 
both ZW and Zl”! regions. However, because of the preferential sampling of the tur¬ 
bulence along the particle trajectories, the drift into Zl+1 regions is much stronger than 
that into Z[“l regions. These results clearly demonstrate that the non-local clustering 
mechanism that dominates St ^ 0(1) particles is fundamentally different to the tradi¬ 
tional centrifuge mechanism, yet the particles continue to drift more strongly into Z^"''! 
regions than Zf“l regions, consistent with the phenomenon of preferential concentration. 

We now use the DNS data to test the key assumptions made in the theoretical analysis 
in First, we made the assumption that for finite St, the Lagrangian autocorrelations 
of S{xP{t),t) and TZ{xP{t),t) are positive, although strictly speaking the analysis only 
requires that they remain positive for max[t — t', t — t"] ^ 0{Tp). Figure shows the 
autocorrelations 

(3.5) 

(3.6) 


^s(t,t') = (‘5(a;P(t),t) : S{xP(t'),t')'j j(s{xP{t),t) : S{xP(t),t) 
^nit,t') = (n{xPit),t) : n{xP{t'),t')) /{n{xP{t),t) : n{xP{t),t)), 


for particles with St = 0(1). The results clearly demonstrate that although 4>5(t, t') and 
are affected by the particle inertia, they remain positive for t — t’/Tp ^ 0(1), 
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Figure 3: Plot of DNS data for ^'5 (solid lines, open symbols) and 'i’-ji (dashed lines, 
filled symbols) for various St. 

confirming the assumption. A second assumption made in our analysis was that the auto¬ 
covariances {Z^{t', t”))z[+-] and {y^{t', t"))z[+,-] remain of the same sign for max[t — t',t — t”] < 0{Tp). 
Figurej^plots the corresponding autocorrelations normalized by their value at t' = t" = t, 
for St = 1. The results indicate that for sufficiently large t — t'^t — t" the autocovariances 
can change sign (see figure]^), reflecting the fact that a particle in a region at time 
t, for example, may have been in a 2 ^ 1+1 region at some earlier time in its path-history. 

However, the results confirm that for max[t — t',t — t”] ^ 0{Tp) the autocovariances do 
not change sign, confirming the assumption in our analysis. 

A third and related assumption made in our analysis is that the timescales associ¬ 
ated with {ZP{t',t"))z[+,-\ and {y^{t',t"))z[+-\ are ^ Oijp) which we estimated to be 
valid based on the known behavior of ts and t-ji . The results in figure show that for 
max[t — t', t — t"] ^ 0{Tp), the autocorrelations remain signihcant (> 0.5) demonstrat¬ 
ing that over the memory timescale of the particle, the strain and rotation fields remain 
significantly correlated along the particle path history. 

These results validate the key assumptions made in our analysis. 


4. Conclusions 

By means of theoretical analysis and tests using DNS data we have demonstrated that 
the non-local clustering mechanism acting in the dissipation range of turbulence with 
finite ts and t-jz is consistent with the phenomenon of preferential concentration. This 
occurs because the particles preferentially sample regions of high-strain and low-rotation 
along their path-history, and this preferential sampling affects the non-local contributions 
to the drift velocity generating the clustering. The St = 0{1) regime is quite distinct 
from the S't <C 1 regime, where preferential concentration is easily understood because 
the mechanism generating the clustering arises from the local preferential sampling of 
the fluid velocity gradient held. 
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Figure 4: Contour plots of DNS data at St = 1 for (a) {ZP{t',t''))z[-]/{2P{t,t))2:[-], 
(b) {ZP{f,t"))zM/{ZP{t,t))zM, (c) {ynt',n)zi-,/{yPit,t))zi-, and (d) 

{ynt'X))zi+}/{ynt,t))zw- 


For a flow with rg —>• 0 and t-ji —>■ 0, such as in the model flow field used in Gustavs- 


son & Mehlig (20111, or in Navier-Stokes turbulence with particles settling under strong 


gravity, the centrifuge mechanism no longer operates and preferential concentration van¬ 
ishes. Nevertheless, inertial particles in such a system cluster because of the non-local, 
path-history symmetry breaking mechanism. 
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